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Abstract 



We compute the free energy density for pure non-Abelian gauge theory at 
high temperature and zero chemical potential. The three-loop result to 0{g^) 
is 
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where T is the temperature, C, is the Riemann zeta function, g^ = g{li)C^ , Jl 
is the MS renormalization scale, g{ji) is the corresponding coupling constant, 
and c^A and Ca are the dimension and Casimir of the adjoint representation. 
We examine the sensitivity of this result to the choice of renormalization scale 
Jl. We also give a result for the free energy of scalar (f)^ theory, correcting a 
result previously given in the literature. 
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I. INTRODUCTION 



The perturbative expansion of the free energy of hot non-Abehan gauge theory is of the 
form 

F ~ T^[co + C2g' + c^g' + (4 In^ + c^)g' + 0(/)] , (1.1) 

where the Cj are numerical coefficients (with some dependence on the choice of renormal- 
ization scale). The leading term is just the free energy of an ideal, ultrarelativistic gas of 
bosons. The first effect of interactions appears at 0{g^) and can be computed from two- 
loop diagrams such as fig. |l|. To compute to higher order requires reorganizing perturbation 
theory to account for Debye screening of electric fields in the plasma and yields terms non- 
analytic in g^ such as 0{g^) and Olg"^ In g). The full 0(5'^) term requires a 3-loop calculation, 
and a full accounting of Debye screening at 3 loops would produce the 0{g^) terms. And 
that's it; perturbation theory is believed incapable of pushing the calculation to any higher 
order. Beginning with four loops, infrared problems associated with magnetic confinement 
appear and non-perturbative 0{g^) contribution to the free energy.[] A complete three-loop 
calculation of the free energy therefore has the special significance that it's the best anyone 
will ever do with perturbation theory. In this paper, our goal is slightly more modest. We 
shall only tackle the 0{g'^) contribution from three loops and leave the 0{g^) contribution 
for another day. 

Another interest of the three-loop calculation is that 0{g'^) is the first order that begins to 
implement the renormalization-scale independence of the free energy. The coupling in ( |1 . 1|) 
is really g{^) where /x is some renormalization scale, and some of the coefficients depend 
on In/i. The leading term that depends on the interaction is order g'^{fi), and by itself 
depends logarithmically on our choice of /i. A change in this term due to a small change in 
renormalization scale. 



^ For a review of this, and also of the previously mentioned reorganization of perturbation theory due to 
Debye screening, see sec. IV of ref. and also ref. 
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FIG. 1. A 2- loop contribution to the free energy. 



g\^,')=g\^)+P,g\^^)\n^ + ... ^ (1.2) 



is compensated by changes in higher-order contributions, first starting at 0{g'^). The 0{g^) 
result should therefore have a flatter dependence on /i than the 0{g'^) result. By checking 
this claim, we can get some idea of the theoretical uncertainties of lower-order calculations 
and perhaps learn some qualitative lessons that will carry over to other thermal quantities. 

The 0{g^) piece of the free energy of non-Abelian gauge theory was previously obtained 
by Kapusta 0, and the 0{g^ liag) piece by Toimela We shall compute an analytic result 
for the full 0{g'^) contribution. In somewhat related work, Coriano and Parwani |^ have 
recently studied high-temperature QED and numerically extracted the 0{g^) contribution, 
and Parwani [^] has also found the 0{g^) piece. (Unlike in non-Abelian gauge theory, the 
perturbation series in QED does not break down after g^.) We shall only study pure gauge 
theory in this paper and do not include any fermions. Fermions will be included in a later 
work. 

In the next section, we warm up to our task by computing the 0{g^) contribution to the 
free energy in pure scalar theory. The result for the basic, three-loop scalar diagram will be 
essential to the later gauge theory calculation, and we shall step through our technique for 
evaluating it analytically. We shall also briefly review the reorganization of the perturbation 
theory to account for the scalar analog of the Debye mass. In section |T|, we turn to non- 
abelian gauge theory and show how many 3-loop diagrams can be reduced to the scalar 
case. We then discuss how to evaluate the exceptions, which are two-particle-reducible 



diagrams. Finally, in section IV we discuss our results and examine the renormalization 




FIG. 2. 1-loop contribution to the scalar thermal mass. 

scale dependence. The details of several calculations needed along the way are relegated to 
appendices. 

Throughout this paper we shall find it convenient to work almost exclusively in the 
Euclidean (imaginary time) formulation of thermal field theory. We shall conventionally 
refer to four-momenta with capital letters K and to their components with lower-case letters: 
K = {ko, k). Unless explicitly noted otherwise, all four-momenta are Euclidean with discrete 
frequencies /cq = '^tttiT. 

II. SCALAR THEORY 
A. Basics 

Consider the theory of a real-scalar field with Euclidean Lagrangian 

>Ce = ^(50)^ + ^^7V, (2.1) 

where we consider temperatures large enough that any zero-temperature mass can be ne- 
glected. The 0{g'^) contribution to the free energy of this theory has been computed nu- 
merically by Frenkel, Saa, and Taylor In this section, we show how to obtain the result 
analytically and also correct an error in the derivation of Frenkel et al. 



^ Note that our is 4! times their and that our e is half of theirs. 
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At high temperature, the scalar picks up a thermal contribution ^g^T"^ to its effective 
mass from the one-loop diagram of fig. It is inefficient to do perturbation theory with zero- 
temperature scalar propagators, which do not account for this effect. We follow refs. 
and rewrite the Lagrangian a^ 

C^ = C, + \g^' - ^^g'T^' , (2.2) 



(2.3) 



where the thermal mass has simply been added in and subtracted out so that nothing is 
changed. Now treat £o as the unperturbed Lagrangian and the last term as a perturbation. 
This reorganization of the perturbative expansion is necessary to get a well-behaved expan- 
sion in g. One can imagine including yet-higher order corrections to the thermal mass in Cq 
above, but this is unnecessary and we shall not do so. 

We regularize the theory by working in = 4 — 2e dimensions with the modified minimal 
subtraction (MS) scheme. This corresponds to doing minimal subtraction (MS) and then 
changing the MS scale to the MS scale Ji by the substitution 

In dimensional regularization, the one-loop thermal mass generated by fig. is 



'n' = \9%j,, (2.5) 
where the integral-summation sign above is shorthand for the Euclidean integration 

t-^-TY:l^ (2.6) 



PO 



and the sum is over po = 27rnT for all integers n. Our reorganized Lagrangian is 

Ce = Co + \{Z, - + j/'Z.g'ct^' - \m^^ , (2.7) 



(90)^ + m^0^ , (2.8) 



^ For a short review in a slightly different context, that also contrasts this resummation scheme with the 
slightly different one we shall use in the next section for gauge theories, see ref. fl. 
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(e) 





(h) 



FIG. 3. Diagrams contributing to the free energy in scalar cj)'^ theory. Dots represent zero-temperature 



counter-terms, and crosses represent the "thermal counter-term" arising from the last term of (2.7) 



where Zi and Z2 are the usual zero-temperature multiphcative renormahzations: 



Z, = l + Oig'), 
3 9' 



1 + 



+ Oig'). 



(2.9) 
(2.10) 



2e (47r)2 

The diagrams contributing to the free energy F through three loops are shown in fig. 
where all propagators represent the reorganized propagators of Co- The sum of these di- 
agrams give —F. All of the diagrams except the last, the basketball diagram, are simple 
because they factorize into one-loop integrals. Diagrams (e-g) are particularly simple be- 
cause they cancel each other at 0{g'^). Diagram (a) represents the contribution to —F of a 
non-interacting gas of bosons of mass m, and its high-temperature expansion is well-known 



^ ln(p2 + ^2 



^ rr4: ^ 2rr2 , 3rr , ^ 

— I ml H ml H - 

90 24 127T 647r2 



- + 21n -4^ + 27, 



m 



+0(mVT2,e) 



(2.11) 



The one-loop integral needed for the remaining diagrams is obtainedPl by differentiation with 



^ For a sketch of an alternative derivation directly in Euclidean space, see for example section III.D of 
ref. §. 

6 



respect to m : 



1 



'p P"^ + m? 



— (1 + eL,)T^ — -mT - 
12^ ^ 47r (47r)2 



m 



i + 2111-4^ + 27, 
L e 47rT 



(2.12) 



We have shown only those terms that can contribute to the free energy at 0{g^), but this 



requires introducing a term of 0(e) that wasn't needed in (|2.11|) . The coefficient of this 
term is less well known, and it is worth taking a moment to focus on the m = case and 
review its simple derivation: 



1 

P P2 



^ J f27r)3-2. ^ f27r)3-2^ 



(27r)3-2^p2 



2/i2^r J 

2v^ \nT^ 
T2 



2e 



P 



n=l 

oo 



(27r)3-2^ (27mT)2+p2 



n 



(27r)3-2^ (27rT)2 + p 

2 / 2 \ ^ 

r(-l + e)a-l + 2e 
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SO that 



2 In 



47rT 



/i 



c(-i; 



2^ + 2. 



(2.13) 



(2.14) 



AnT ' C(-l) 

With these tools, all of the diagrams but the last are straightforward. The individual 
contributions of each diagram are summarized in appendix 0. Frenkel et al. [0 did not 
properly account for the term of ( p^.l2[ ) when evaluating diagrams (a-c)J^ 



B. The basketball diagram 

The last diagram of fig. ^ is more difficult. One simplification occurs because the diagram 
remains infrared convergent when the mass m in the propagators is set to zero. This means 



^ More specifically, because their definition of mn? = ^9^T^ differs from ( ^.5[ ) at 0(e), their thermal 
counterterm does not exactly cancel the one-loop diagram of fig. ^ and they should have an extra term in 
their eq. (14). If one instead uses our definition ( ^.5[ ), then their eq. (14) is correct but there should be an 
extra term in the one- loop pressure in their eq. (11). 



7 



that the presence of m has only a sub-leading effect on the diagram, since m itself is 0{g). 
We can ignore m here if we are interested in the free energy only to 0{g'^), and the basketball 
diagram is then proportional to 

= tpQK pWkW+Q + ky • 

Our attack on this integral starts with the observation that the basketball diagram of 
fig. |^(h) requires three independent four-momentum integrations if evaluated in momentum 
space but only one four-space integration if instead evaluated in configuration space. This 
suggests the diagram may be more tractable in configuration space. Indeed it would be if 
the diagram were ultraviolet (UV) convergent and we could set e to zero. The configuration 
space propagator in four dimensions, with period 1/T in Euclidean time, is the relatively 
simple function 

A(r,r) = ^, sinh(2vrrT) ^2.16) 

ATtr [cosh(27rrT) — cos(27rrT)J 

but in 4— 2e dimensions it is a nightmare. We therefore need to first subtract out the UV 
divergent pieces, and evaluate them separately, so that we can then evaluate the remainder 
in four dimensions. We found, however, that making these subtractions is more convenient 
in momentum space than in configuration space. As a result, our derivation mixes the use 
of momentum and configuration space. First, we shall always treat the Euclidean time 
direction in frequency space. In configuration space for the remaining, spatial dimensions, 
the propagator 1 /P^ then has a very simple form in exactly four dimensions: 

A(j>o,r) = ^. (2.17) 



Our approach will be to start with the the momentum space form ( p.l5| ) of the basketball in- 



tegral, convert it step by step into a configuration space form, and make needed subtractions 
as we go along. 
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1. A careless derivation 



To simplify the presentation, let's first forget about the UV subtractions and run through 
the derivation pretending that it makes sense in exactly four dimensions despite the UV di- 
vergences. We'll later step through it a second time, handling the divergences more carefully. 
We start by noting that, by a shift of variables, the momentum integral ( |2.15| ) can be written 
in the form 



'ball 



^^[n(P)]2, (2.18) 
where 

Now let's evaluate n(P) using configuration space and ( |2.17| ): 
n(P) =TY,(d'r e^P"^'A(go, r) A(po + go, r) 

90 

T r 1 

= / d^j. ^iV-r ^-\qo\r ^-\po+qo\r 

(47r)2^i r2 

/ ^^e'^^'lcothf + |po|)e"IP°l' , (2.20) 



(47r)2 

where we have introduced the dimensionless variables 

f = 2'KTr , po = Po/'^T^T . (2.21) 
Plugging this into (|2.18|) , the p integral becomes trivial, producing 



^baii = -J^^Y. J ^^^^ ^(cothf + IpoD'e 



-j^2-2\po\r 



Po 

2^4 »oo 



df (coth^ f - coth f dr + coth f . (2.22) 



327r2 

As we shall discuss later, integrals like this can be performed analytically when they are 
convergent. The present result doesn't make any sense, however, because the UV behavior 
(f 0) of the integrand makes the integral divergent. We shall now repeat the above 
derivation while making necessary subtractions as we go along. 



2. Subtraction of UV divergences 



Let's start with the expression (|2.20|) for n(P). This integral is logarithmically divergent 
in the ultraviolet. As usual with one- loop integrals at finite temperature, however, it can be 
made finite simply by subtracting out the zero-temperature contribution. So we write 

n(P) = n(°) (P) + n(^) (P) , (2.23) 

where n(°)(P) is the zero-temperature result 



In four dimensions, n*-"^-* (P) can be obtained from ( |2.2U| ) by subtracting out its T — limit 
(with P fixed): 

n(^)(P) = / d^r ^e'^-'' ( cothf - -\ e-|P«l'- + 0(e) . (2.25) 
[An)^ J \ r J 

We now split our computation of the basketball diagram into 

/ball = ^[n(^)]^ + 2^n(^)nw + ^[nW]^ . (2.26) 

Though n^"^) is finite in four dimensions, the first term above is not because 11*^-^) ~ 1/P^ 
as P— s>oo; the first term therefore has a logarithmic UV divergence. The large P {i.e. 
P :^ T) behavior of H^^) is easy to extract by staring at the definition (|2.19| ) of 11. The 



dominant contribution comes from routing the large momentum P solely through one of the 
two propagators and then integrating over the relatively small momentum Q ^ T in the 
other propagator: 

A more rigorous derivation may be found in appendix |B[ This limit is not restricted to four 
dimensions, so we are now in a position to subtract out the UV divergence in our integral: 
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The first integral is now UV finite, and so we might hope to evaluate it in exactly four 
dimensions. However, it is not also infrared finite if we evaluate the Po = term of the 
frequency sum; for po = 0, the subtraction we made diverges linearly with pin the infrared. 
We shall therefore treat the po = mode separately and put primes on integrals, as we have 
above, to denote that this mode is excluded: 



2e 



P 



\3-2e 



(2.29) 



Po^o^ (2-) 

We can now evaluate the first term of ( p.28|) in exactly four dimensions. 

The leading large P behavior ( p.27|) of 11*^^^ is related to the leading small r behavior of 
the integrand in ( p.25|) and is given by 

T /• 1 



(2.30) 



Following the same steps as in the careless derivation of the previous section, we then obtain 



3^ 



coth r — - 
r 



l\2 /f\2 



3) 



rp4 



327r2 Jo 



drr~ 



coth r — 



g-2|po|r- ^ 0(e) 

(cothf- 1) + 0(e). (2.31) 



This integral is both IR and UV convergent and can be evaluated using the techniques of 
appendix to give 



^^'|[n(^)(p)p-(^ 



-i - 



1 (T 



2\ 2 



C'(-3) , C'(-l) 28 



(47r)2 \^12; 

The last term in ( |2.28| ) is easily evaluated in 4— 2e dimensions using ( ^.13| ) and 



+ 0(e). (2.32) 



1 , u 
- + 21n/- + 27. 



fp P4 (47^)2 

which may be obtained in a manner similar to ( p. 13 ). 



0(e) 



(2.33) 
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Completing the derivation of the [11*^-^)]^ contribution to /baii now just requires adding 
in the contribution of the Po = mode, which is UV convergent and does not require any 
subtractions: 



T 



327r2 JO 



CO / \' 

^^-'i I cothf — - ) + 0(e) 
r, 



dr r' 



(i^(S) |-"«C'(-2)|-H0W, 



(2.34) 



where we have again used the techniques of appendix |C| to do the integral. Putting together 
( CT , iH, H, m, Ipl) then gives 
1 



'ji2 \ 2 r 



^ + 241n^-16^ + 40^ + ^ 
e A-kT C(-3) C(-l) 15, 



+ 0(e) 



(2.35) 



(47r)2 \^12 

Now that we've covered the basic ideas of our technique, we'll leave the evaluation of the 
remaining two terms in ( p.26| ) to appendix 0. The final result for the basketball integral 
( p7[5| ) is 

1 



'ball 



(47r)2 1,12^ 



^ + 361nJ^-12^ + 48^ + l^ 
e AtcT C(-3) C(-l) 5 



OU). 



(2.36) 



This agrees with the numerical result of ref. We should mention that our analytic result 
can also be obtained from the integrals generated by a real-time analysis, such as in ref. [0 , 
and we show how to do this in appendix We have found it simpler to stick to Euclidean 
space, however, to evaluate diagrams involving double poles l/P"^ which will appear later in 
gauge theories. 

C. The result 



Putting together all the diagrams, which are independently tabulated in appendix 
one finds that the free energy in 0^ theory at high temperature is 



rp4 



9 1 10 8 V47ry V6 V47ry 



3, lC'(-3) IC'(-I) 1 59 



AttT 4C(-3) 2C(-1) 



120 



+0(r 



(2.37) 
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FIG. 4. The setting sun diagram. 

Before leaving scalar theory, we should mention one other basic scalar integral which 
appears in the literature [^,0 and which will be needed when we analyze gauge theories. It 
is the integral corresponding to the setting sun diagram of fig. 

Isnn{'mi,m2,'m-i) = ^ 2..„2 , 2\\tTD , r^\2 , 2T ' (2.38) 

TpQ [P^ + mf}{Q^ + m^) [{P + QY + mil 

evaluated to leading order in the masses. It has previously only been evaluated numerically 
0], but the same techniques we applied to the basketball diagram can be used to obtain an 
analytic result. We give the derivation in appendix with the result that 



(47r) 



1 + f 



4e \mi + 7712 + J 2 



0(771, e). (2.39) 



III. NONABELIAN GAUGE THEORY 

We now turn to pure non-Abelian gauge theory, given by the Lagrangian 

Ce = \ {d,A: - d.Al + gr'^AlJ^)^ + (gauge fixing) . (3.1) 

We shall work exclusively in Feynman gauge. (It would be nice to explicitly verify that our 
results are independent of gauge choice, but we have not done so.) Let dp^ and Ca be the 
dimension and Casimir of the adjoint representation, with Ca given by 

For SU(A^), they are 
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FIG. 5. The one- loop gluon self energy. 
dA = N^-l, Ca = N. (3.3) 

It is also convenient to define the effective coupling of the adjoint representation by 

gl^g'CA- (3.4) 

As before, we shall regulate the theory with dimensional regularization in the MS scheme. 

As in the scalar case, one-loop effects induce a thermal mass contribution. This mass is 
given by the one-loop self-energy 11^,^ at zero momentum, and in Euclidean space a mass is 
generated for Aq but not for This mass M is the Debye screening mass for static electric 
fields and may be evaluated from the diagrams of fig. |^ as 



= Uf^iO) = gl{d - 2ff^ ^^6'^' . (3.5) 



In four dimensions, is simply g\T'^/'i. In the gauge theory calculation, we find it 
calculationally convenient to use a slightly different reorganization of perturbation theory 
than we did for in the scalar case. The success of the reorganization only depends on the 
behavior of the propagator in the infrared (po=0, p<^T), where the mass cannot be treated 
as a perturbation. We follow ref. and only introduce the mass for the po = mode. That 
is, we rewrite our Lagrangian density, in frequency space, as 

£e = (-Ce + \M^AIA15,^) - \M^AIA15,^ , (3.6) 



^ Throughout this article, 11^^(0) will refer to the Euclidean limit (po=0,p-^0) and never to the limit 
(j)Q-^Q,p — 0) which may be achieved by analytic continuation and which gives the mass gap for propagating 
plasma waves. 
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where Sp^ is shorthand for the the Kronecker delta function Sp^^. Then we absorb the first 
term into our unperturbed Lagrangian Cq and treat the second term as a perturbation. [| 
Since the necessity of resummation is an infrared phenomena, associated with the mass 
scale gT, it is worth noting that the prescription ( p.6| ) can be naturally expressed in the 
language of decoupling. First imagine integrating out all the physics associated with scales 
^ T. In particular, integrating out all of the Po ^ modes in Euclidean space generates an 
effective three dimensional theory of the remaining po = modes. This effective theory will 
have the thermal mass for Aq and other interactions induced by the heavy modes, which 
can be computed to any desired order in perturbation theory. Only then does one finally 
integrate out the po = modes after deciding on a sensible partition of the effective three- 
dimensional Lagrangian into an unperturbed piece, containing the thermal mass terms, and 
a perturbative piece. Rather than carry out the reduction to this effective theory explicitly, 
however, we find it simplest to just introduce the reorganization ( |3.6| ). We refer the reader 
to sections III.D and VI of ref. |p for details of how to implement this form of reorganization 
on two-loop graphs.^ 

The diagrams that contribute to the free energy are shown in fig. ^. The diagrams 
involving only one-loop integrations are trivial, and the resummation and the two-loop 
graphs can be handled by the methods of ref. [§. Let's therefore focus on the three-loop 
diagrams. The first potential problem is that some of the individual diagrams, such as 
fig. 1^, are infrared divergent because of the masslessness of the A propagator. However, the 
particular combination we have shown in fig. ^(1) is well-behaved in the infrared since the 
shaded blobs. 



^ This reorganization only helps in the evaluation of static quantities such as the free energy. To evaluate 
time-dependent correlations in real time, one would need the resummation scheme of Braaten and Pisarski 

^ The reader should beware that many of the specific formulas of ref. are particular to Landau gauge, 
whereas in the present work we are working in Fcynman gauge. 
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(a) 



(f) 



(b) 




(g) 




(k) 



(c) 



(h) 



(1) 




(d) 



(i) 



c 



V 



(m) 




(e) 




(j) 



FIG. 6. Diagrams contributing to the free energy in gauge theory. The crosses are the "thermal 
counter-terms" arising from the last term of ( |3.6| ). We have not explicitly shown any zero-temperature 
counter-terms, and each diagram should be multiplied by the appropriate multiplicative renormalizations 
for vertices and propagators. 
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FIG. 7. A diagram with infrared problems. 

An^, = n^,(P) - n^.(o)5p^ , (3.7) 

are 0{pT) for {po=0,p—^0). We can then also drop the Debye mass M in evaluating the 
three-loop graphs since, as in the scalar case, the corrections to the free energy will be 
beyond 0{g'^). 

All of the three loop diagrams except (1) can be reduced to the scalar basketball integral 
of ( p.l5 ). For example, diagram (i) is equal to 



_ _ _1 , P-{Q-K) {P-K).Q 

This may be reduced by (1) expanding numerator factors in terms of denominator factors 
to cancel factors between numerator and denominator, such as 

P-{Q-K) = \[{K- Pf -K^-{P- Qf + g2] . (39) 

(2) performing appropriate changes of variables to collect similar terms; and (3) using the 
identity 



(3.10) 



rp{P + Qf{P + Kf 2 Tp{P + Qf{P + Kf' 

which follows by averaging the left-hand side with itself, after applying the change of vari- 
ables P —P—Q—K. Appendix |^ steps through this reduction for the example ( |3.8| ), and 
the reductions of diagrams (g-k,m) are all tabulated in appendix ^ 

Unfortunately, diagram (1) cannot be reduced to the scalar basketball. If one tries the 
above tricks, one finds a term of the form 
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TpqkP^Q^K\P + QY{P + KY ^ ■ ' 

for which the tricks fail to remove the numerator factor. So we have a new basic integral 
that we must evaluate, like the basketball integral of scalar theory. We have found it more 
tractable, however, to apply our integration method directly to the original diagram (1) 
because the orthogonality of the one- loop self-energy Ii^j,{P) to P^ leads to useful algebraic 
simplifications. Diagram (1) is proportional to 

dAQihcA = -^tr[An^,(P)]2 . (3.12) 

The evaluation of Jqcd is somewhat similar to that of the basketball integral Jbaii and is 
presented in appendix 0. We should mention, however, that the derivation is more compli- 
cated and involves a miraculous cancelation between two complicated integrals that we don't 
know how to calculate individually. The appearance and cancelation of such complications 
suggests that we may still be missing the most elegant method for making these calculations. 

All the results for individual graphs are collected in appendix with the final result 
that 

F = ciAT^^|-i+ I'A^-^ f^V -48 m\J^] 
9 1 5 V47ry ^W) W) l27rv^/ 



+ 



9aV 



^1 38 C(-3) 148 C(-1) , 64 

3 ^ AttT 3 C(-3) 3 C(-l) 5 



AttJ 

For those who prefer C, functions with positive arguments. 



+0{A)\- (3.13) 



ft:4^.„(2.)+,.-i:i-fii±4. nodd. (3.14) 

C(-n) C[^ + n) 



IV. DISCUSSION 



Evaluated numerically, our result ( p.l3|) is 



F = -dA- 



45 



1 -0.31250 ( — ) +0.7216 



TT 



71 J 



+ 



9A 



TT 



0.93750 In — - 0.14323 In ^ + 0.74582 

TT T 



(4.1) 
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thru 
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thru ^ 

thru ^ excluding g 3 
thru § log |l 



2 4 6 

loao ( P/T) 

FIG. 8. The dependence of the free energy F on the choice of renormahzation scale jl for pure gauge 
QCD with as{T) ~ 0.1. The free energy is norniahzed in units of the ideal gas result — (iA7r^T''/45. The 
thick solid, dashed, and dotted lines are the results for F including terms through g"^, g^, and g^ respectively. 
The light solid curve is the g"^ result plus the g^ ln(/2/r) term required by renormahzation group invariance. 
The light dashed curve is the g"^ result minus the g^ term. 

We have chosen the expansion parameter gA/i^ simply because it makes all the coefficients 
0(1). 

Now we can ask whether perturbation theory is behaving well for physically-realized 
values of the couplings. In particular, we can investigate (1) the size of corrections from 
different orders for a fixed choice of renormahzation scale, such as /i = T, and (2) whether 
higher-order results are less sensitive to the choice of the renormahzation scale fl than lower- 
order results. This information is summarized in fig. ^for pure gauge QCD with as(T) = 0.1, 
which for real QCD would correspond to a temperature around the electroweak scale. We 
have used the two- loop renormahzation group to compute g{fi): 

^^^-/3oln^ + |iln(l-/3o<?I(T)ln^) , (4.2) 

where 
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'3(47r) 



68 

'3(47r) 



(4.3) 
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FIG. 9. The same as fig. || but for as{T) = 0.02 [which is identical to SU(2) gauge theory with 
aw(r)«l/33]. 

At /2 = T, the terms of ( [4.1[ ) for as (T) =0.1 are 

F = -dA——[l - 0.12 + 0.17 + (-0.07 + + 0.11) + 0{gl)] . (4.4) 
45 

The behavior of the perturbative expansion doesn't look particularly good, though a partial 
cancelation between the g'^lng and g'^ terms makes the total 0(5'^) contribution relatively 
small. Alternatively, examine the sensitivity of the result to the choice of renormalization 
scale by examining the slopes of the curves in fig. |^ at /i = T. Contrary to one's expectation 
for a well-behaved perturbative expansion, the 0{g^) result is more sensitive to ft than the 
Olg"^) or 0{g^) results. 

The 0{g'^) result has to be less sensitive to /2 if (? is sufficiently small. Fig. ^ shows the 
dependence for as{T) = 0.02. This is equivalent to a system of interest — pure electroweak 
theory at the electroweak scale, with ^ 1/33. Yet still the 0((7'^) result is no less 
sensitive than the 0{g^) result. Fig. ^ shows as{T) = 0.001, where we finally see the 
expected behavior. (In six-fiavor QCD, this as would correspond to a temperature of 10^^^ 
GeV.) 

The source of the sensitivity problem can be found by remembering that the g^ term 
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FIG. 10. The same as fig. | but for UsiT) = 0.001. 

probes different physics than the term: the term is produced by particles with hard 
thermal momenta of order T, while the g^ term arises from the interactions of particles 
with softer momentum of order gT. Since the g^ term is the leading-order contribution of 
the physics of scale gT, it perhaps should be treated independently of the g^ term when 
discussing whether the perturbation expansion is well-behaved. The g^ terms reduce the 
sensitivity of the g^ term to /i, and the g^ contribution will be needed to reduce the sensitivity 
of the g^ term. The light dashed line in figs. ^|-^ show the result of the free energy through 
0{g^) if the g^ term is artificially excluded. The sensitivity to /i, compared to the 0{g'^) 
result, is indeed much better than before. 

In order to put the 0{g^) term back in, we have tried adding the g^lnfi term that's 
determined by renormalization group invariance. The light solid lines in figs. represent 
our result ( |3.13|) with the addition 



5 r 



11 In 



(4.5) 



9 [ ^/3 \AnJ 

The results are much better behaved than the 0{g^) results discussed earlier. Of course, the 
constant under the log at 0{g^) is unknown and will change the curves somewhat. Our under- 
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standing of whether perturbative results are indeed well-behaved in high-temperature QCD, 
for realistic coupling constants, would therefore benefit by a true calculation to 0{g^) — the 
last order accessible to perturbation theory. 

This work was supported by the U.S. Department of Energy, grant DE-FG06-91ER40614. 
We thank Larry Yaffe, and most especially Lowell Brown, for their repeated suggestions 
that we try analyzing our basic diagrams in configuration space using the simple form of 
the propagators in three spatial dimensions. We also thank Rajesh Parwani and Claudio 
Coriano for useful discussions. 
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APPENDIX A: RESULTS FOR INDIVIDUAL GRAPHS 



1. Scalar theory 



The diagrams of fig. § are given by 



-f^''F^ = -l9'(iT,T^)\ (Alb) 



+ + = 0(/), (Ale) 
V^i^' = ^Aau + 0(/). (Alh) 



2. Gauge theory 



Writing F = fi "^"^dxT and ignoring terms of 0(e), the diagrams of fig. |^ are given by 



2T 127r 
lnP2, 



- — (M^ + + Ml)MT , 



2 ry2 



9iz, 

Si 

52 



p2 



+ 5i + 52 



1 MfMT 



resum ~l~ o 2 



■4 (47r)2 



p2 ' - 



+ —M^MT . 



57r 



^rf-|)^i4all + 0(/) 



(A2a) 
(A2b) 
(A2c) 

(A2d) 

(A2e) 



(A2f) 
(A2g) 
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-^' = -^^?i4all + 0(/) 

Id 



32 



^i/ball + 0{i 



-^j = -rf(rf-l)(7iJball + 0(/) 

27 
lo 



:^7l/ball + 0((7^ 



(A2h) 
(A2i) 
(A2j) 
(A2k) 
(A21) 
(A2m) 



The multiplicative renormalization constant used for the coupling is given by 

11 9l 



5'ba 



^99 



9 



(A3) 



6 (47r)2e 

Wave function renormalization constants are unnecessary because they cancel between ver- 
tices and propagators. [To this end, the most convenient choice of M is 



(A4) 



which differs from ( p.5| ) by the introduction of Z| and the photon wave function renor- 



malization Z^. At 0{g ), however, this is not an issue — the factor of ZgZ\ in M can be 
ignored for all of our diagrams.] , M|, and M| denote the three pieces of originating, 
respectively, from the three diagrams in fig. |^. We will not give explicit formulas for these 
pieces because they explicitly cancel between diagrams (c-f). 

The integrals needed above are given by ( pTj |2]T|, |]36|, |ATO| , 0, The 
integral /resum is defined by 

















[P + QY Q'\ 










+ 0{g') 




Q? 





(A5) 
(A6) 



and will be discussed below. 

The effect of the thermal mass term in the po=0 gauge propagator appears fairly simply in 
the last terms of ([A2| a, d) and in (|A2| c) . The case of diagram (e) is a little more complicated. 
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The first term of (|A^e) is tlie result wlien M is ignored. 5i is tlie correction to tliis result 
for the contribution to the diagram where exactly one of the three propagators has po = 0, 
and 82 is the correction for the contribution where all three have po = 0. So, for instance, 



5i 



AT {p^ + M^)Q^{p + qy 



[{P - Q)^5ou + 2Qo5^, - (2P + Q)J^of ~ {M-^0} 



(A7) 



where the factor in brackets is the triple gauge vertex. Using the reduction tricks described 
after (|3.8|), 61 may be reduced to 



(5i = -M^ 



+ {d-l) 



+ 



1 



<il 



Q\P + Qf 



(A8) 



(p2 + JVf2) Q2 ' 

which may be recast in the form shown in (|A2| e). ^2, and the mass effects in (|A^), are 
calculated similarly. 

-^resum is casily evaluated by taking the form (|A6|) and scaling all three- momenta by j^ol' 



go 



X 



1 



+ 0{g'). (A9) 



;i + g2)(l + |p'+g-]2) (l + g2)2_ 

Recognizing the go sum as giving a ^-function, and that the integrals are finite, it is easy to 
now take the e— i>0 limit: 



2M^TX(0) 



dP'p d?q 1 
(27r)3(27r)3^ 



1 



;i+g2)(l + |p+g]2) (l+g2)2 



0{g\e) 



8 (47r)^ 



+ 0{g\t). 



(AlO) 



APPENDIX B: LARGE P BEHAVIOR OF SCALAR n(P) 



In this appendix, we shall derive the large momentum behavior of n(P). Recalling the 
definition ( pAOD , 
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one may sum over the bosonic frequency modes by the usual contour integral trick.0 Sub- 
tracting out the zero-temperature piece gives the finite-temperature part of n(-P): 



n(^)(p) = 



2e^ ^^^^ 



3-2e 



1 



+ 



{q -ipo)^ - \q + p\^ {q + mY-\q-p 



where we have defined the thermal bosonic factor 

n{q) = 



(B2) 



(B3) 



ei/T _ 1 

The exponential fall-off of n{q) for large q ensures that only g ^ T is important. For P ^ T, 
we can then expand the denominators in (|B2| ): 

^3-2e^ n(g) 



{2n] 



3-2e 



2 8{ipoq + p- qf 

p2 ~^ p6 



(27r)3-2^ g 



+ o{Typ^ 



3-2e 



(B4) 



The momentum integrals we need are of the form 



. (27r)3-2^ 

2 / 47r/i^ \ " roo 



(47r)3/2r (I - e) V 
/47r/i2 Y r(3-2e+a)C(3-2e+a) 
V^j 47r3/2r (I - e) 



e'J — 1 



(B5) 



and our large P expansion is 



(B6) 



We note here that 



(B7) 



and so the limit ( p.27| ) in the main text is the same as the leading piece of ( P6|) above. 



See ref. H for a review. 
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APPENDIX C: INTEGRALS OF HYPERBOLIC FUNCTIONS 



In this appendix, we discuss how to evaluate convergent integrals of the form 

1=1 dx\y^ Cm^x"" coth'^ X + > ^ rf™,x'"e-"'"^ 1 , (CI) 



Hdx (^mnX"" COth'^ X + J] rf^x'^e""'"^ j 

\m,n m / 



where a^, c^n, and dm are constants. We shall evaluate such integrals term by term. Since 
the individual terms will in general be divergent, we first regulate in a manner similar to 
dimensional regularization by rewriting I as 



/•oo I 

I = \im dxiy^ Cmnx""^^ coth" x + y^dn 

^ \m,n m 



(C2) 



5 will be used to independently regulate both the x^Q and x— >oo pieces of the integrals, 
again similar to dimensional regularization. Now we can compute three basic regulated 
integrals: 



(C3) 



oo /"OO 

dx coth X = dxx^ 
Jo 



-2kx 



2-^r(z + i)c(^ + i) 



dxx^e 



z„—ax „—l—z 



a-'-'Til + z) 



The rest of the integrals we need are obtainable recursively by 

1 



oo roo 

dx x^ coth" X = dxx 
Jo 

oo 

dx 





n 



-d^ (coth"-^ x) + coth" 



-2 



X 



n — 1 



x""'^ coth""^ X + x"^ coth"~^ X 



(C4) 



(C5) 



(C6) 
(C7) 



After assembling the individual terms of a particular integral (|C^ ), it is straightforward to 
expand in 5 and take the limit 5— >0. 



APPENDIX D: COMPLETION OF THE CALCULATION OF /ball 

In section [II B| , the basket ball diagram was split into three terms, 
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'ball 



(Dl) 



and the first term was evaluated. Here we shall evaluate the two remaining terms. 

We first calculate the second term in ( pi| ). To apply the calculational method of section 
11 B|, we must first subtract the ultraviolet divergences. 

n(o) 

is given by 



n(°) = A 



p2 



(D2) 



where 



A 



^ ^- + 2-7,) +0(e). 



(47r) 



(D3) 



Because Il^'^\P) ~ at large momentum P, the second term of (pi|) is quadratically 



divergent and so requires two subtractions. Using our result ( P6|) for the large momentum 
expansion of we rewrite this term as 



where 



^^n(^)nw = 4 + /b + /e 



n(o)(P) 



(47r)2e 



[n(-)(P)-ngJ(P) 



(D4) 



(D5) 



(47r)2e 



KnP)-ng?(P) 



(D6) 



and 



^^n(°)(P)ngJ(p) 



ng;(P) ^ 2 J_,- + (1 - ^ 



pI 



(D7) 



(D8) 



Ja is ultraviolet and infrared finite, and so it can be evaluated in c? = 4. Note that we have 
used one less subtraction in ( p8| ) for the Po = mode. 
Using the integral representation (|2.25| ) for 11*^^^ gives 
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T 



(47r) 



In 



+ 2-7, 



(Tr—e 



ip-r 



1 

cothf---- + (l-5pJ- 



■|pok 



+0(e) , 



(D9) 



where we have expressed the large P behavior ( |D8| ) of 11*^^^ (P) in terms of three-dimensional 
coordinate space integrals. In fact, it is easy to see that the P— s>oo behavior of Yi^'^\P) 
in ( p^.25D simply corresponds to the f^O behavior of cothf — 1/f, and the subtractions 
in (|D9|) simply reflect the small f expansion 



coth ^=- + -- 7^ + - 
f 3 45 



(DIO) 



The p integral in (P9|) is trivial for the terms that don't involve InP^: it just gives 5(r), 
which in turn gives zero. The p integral involving the logarithm can be evaluated by first 
writing 



d?p 
(2^ 



In 



d 



p^ + Pq da 



dp 



p^ sinpr / 47r/i^ 



27r2 pr + pi 

d 



a=0 



Am'^r da 



(Dll) 



a=0 



Deforming the contour to wrap around the cut in the upper half of the complex plane gives 



d^p 
(2^ 



In 



Att^^ 1 d 

p"^ + Pq da 



sin(7rQ;) / dqqe 



-qr 



1 



27rr \ r 



bo I 

Po\r 



a=0 



(D12) 



Inserting this result into ( p9D and carrying out the po sum yields 
T2 1 



dr \ ( , 1 r r \ ( r d \ , 

, xo , —< cothf \ 1 (cothf- 1) 

(47r)2 2 7o f3\V r 3 45/ V 2 df / ^ ' 



/ , _ 1 f 
+ cothr — - — 
V r 3 



0( 



(D13) 



The integral can be computed by the method of Appendix B to give 
1 



'2^2 \ 2 



(47r) 



12 



8 CX-3) CX-1) 12 _46 
5 C(-3) C(-l) 5 15 



+ 0(e). 



(D14) 
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We now study Jb defined by (p^). Though the P integration converges, we still need 
to evaluate Jb in (i=4— 2e dimensions because of the overall factor of 1/e. Writing n^^^ as 
n - n(°) gives 



1 1 



(47r)2 e Tp 

y4 ^ ' 



nW(P)-ng;}(p) 



(47r)2 e [ 
where we have defined 



^[5o(2)-rf5i(3)] 



S'„(a) = T 



2a-4-2nv/ "yo" 



2 \ ^ ^27r)^+2""2" 



\nT^ J 47r3/2r(a) 
Taking the limit e — 0, 



r (-|+a+e) C(-3-2n+2a+2e) . 



1 /T 



2\ 2 



8 C(-3) I 12 _46 

5 C(-3) C(-l) 5 15 



(47r)2 \12^ 

(We do not have an explanation for the fact that Ja = /b- 
Using the explicit expression ( p2| ) for n*^°\ Jc is 



+ 0(e) 



2A 



47r/i2y f ^ o /I , , 4Ji 

2^2 
2\ 2 



^2 , ^J-i5o(l+e) + ^[5o(2+e)-rf5i(3+e)] 



T 



5e ' 5 "^AttT 5 C(-3) C(-l) 5 + ^ 



(47r)2 \^12 

Adding the results (|D14|) , (|D18|) , and ( P19|) produces 

:n(^)nw 



1 /r2\'^ 



4 24, /i 4C'(-3) X'(-l) 103 

\ In + --^-^^ + 4-^-^^ H 

5e 5 47rT 5 C(-3) C(-l) 15 



(4vr)2 Vl2y 

It is straightforward to also evaluate the third term in f|Dl| 



[n(o) 



47r/i2 

2^2 



2e 



1 fT^V 



So{2e) 



2 12^^ fc 12 C(-3) 24' 

5e 5 AttT 5 C(-3) 5 



(47r)2 \^12^ 

Assembling ( |2.35| ), ( p2CI| ), and ( |D2lD gives our final result ( ^^.361 ). 



+ 0(e) 



(D15) 



(D16) 
(D17) 



(D18) 



+ 0(e). (D19) 



0(e). (D20) 



(D21) 
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APPENDIX E: REAL-TIME CALCULATION OF Ibaii 



In ref. [j^, /baii is expressed as 
1 



'ball 



{Any [12^ 



^ + 181n^ + 187. + 6-36g| 



+ NT^ + 0(e) , 



(El) 



where is given by 

rf^P d^Q d^K 



N = Ar 



{2Txf (27r)3 (27r)3 



327r6 



dp dq dk n{p) n{q) n{k) {p+q+k)\n{p+q+k) — (p+q—k) In \p+q—k\ 



-{q + k — p) hi\q + k — p\ — {k + p — q) hi\k + p — q\ 



(E2) 



and was evaluated numerically to get 



N 



14.17 

327r6 



(E3) 



In (|E2|) , unlike the rest of this paper, P refers to Minkowski rather than Euclidean four- 
momentum, with metric = —Pq + p'^. V denotes that the integrals are to be performed 
with the principal value prescription. The n{p) are the usual Bose factors but in units where 

r=i: 



n{p) 



1 



(E4) 



The second equality in ( p^ ) is obtained by doing the trivial Pq, qo, and /cq integrals and then 
doing the angular integrals. Their final result was then obtained by numeric integration. 
We shall show how to obtain the same result analytically starting from (p^). 

We start by making use of the peculiar identity that {P + Q + can be replaced by 
4|p + g + fcp in the first line of ( [E2| ) for any function n{p) of p = \p\ to give 

^ f d^P d^q d^k n{p) n{q) n{k) 1 

i (27r)3 (27r)3 (27r)3 2p 2q 2k + + k\2 ' ^ > 

This identity can be proved by brute force by doing the angular integrals in (^) by the 
same steps ref. used for (^) and verifying that the result is the same as (^) : 
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N 



d^p f d^q 



(27r)3y (27r)3yo 



'^dk^'^nik) 



2p 2q 



In 



\p + q\ + k 

\p + q\ \\p+ q\- k\ 
\p + q\ + k 



1 r<x> rp+q 

— / dpdqdkn{p)n{q)n{k) d\p + q\\n 
d27r° JO Ap-iI \\P + (l\~k\ 

dp dq dk n{p) n{q) n{k) [p+q+k) \n{p+q+k) — {p + q — k)\n.\p + q — k\ 



327r6 



+ (\p — q\ — A;) In I |p — g| — k\ — {\p—q\+k) ln(|p— 

Sadly, we do not have a more elegant derivation. 

Now take (|E5|) and convert it to coordinate space using the Fourier transform 



(E6) 



d^p n{p) 
(27r)3 p 



Eq. ([E5[) then becomes 



JO 



d r 



dp 



smpr 
eP - 1 



47rr 



coth(7rr) — 



vrr 



(E7) 



1 1 



4(47r; 
1 1 



■ 1 






^coth nr 


Airr 




POO 




do 


dx — (c 




x^ \ 



1 

vrr 



4 (47r)2 
14.1723 
327r6 



47rr 

1x3 

X 



l C(-3) I C(-I) 7 
3 C(-3) 3 C(-l) 45 



(E8) 



Eqs. ( |E1| ) and ( p^ ) give the same result for the basketball integral as (|2.36|) , which was our 
result using the Euclidean formalism. 



APPENDIX F: DERIVATION OF L, 



Let us evaluate the integral for the setting sun diagram defined by (|2.38|) in a similar way 
as the basketball integral. As usual, we work in the limit that masses are all much smaller 
than T, and we shall denote their order of magnitude simply as 0(m). Only the leading 
term in the m/T expansion will be calculated. To this order, the masses can be taken to 
be zero except in the ^0=90=0 contribution, where the mass cuts off a logarithmic infrared 
divergence. But it is convenient to keep only one mass non-zero in the Po=Q'o=0 contribu- 
tion and to set m2=m3=0. The discrepancy introduced by doing so is easily computed in 
coordinate space to be |^ 
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(Pp d?q 1 
(27r)3(27r)3p2 + ^2 



(g2 + ml) [{p + g ^ 



ml 



(47rr)'^ 



-m\r 



'(m2+m->,)r 



+ 0(m,e) 



+ 0(m,e) 



In 



mi 



(47r)2 mi + m2 + ms 



+ 0(m,e) . 



(Fl) 



1. Quick derivation using the contour trick 



We now need to compute 

-^si 1 n ~ 



+ 



1 — 5„ 5a 

^" ^" +0(m,e) 



We start with the purely three-dimensional contribution O 



(F2) 



^ (P2 + r/if) (g2 + ml) [(P + Q)2 + mi] (47r)2 



T2 ri , 

— + In — 

4e mi + m2 + ma 2 



1 



+ 0(e) 



(F3) 



which follows from 



4 



2er 



■^2 / 1 



(p2 + mf)pi+2^ (47r)i-T(l-2e) 
T2 



f— In i 
(47r)2 [4^ ^ ^ 2 



+ 0(e) 



(F4) 



and (pil) . For the second term of (|F2|) , note that if dimensional regularization is used to 
regulate the infrared as well as the ultraviolet then 



4 







r P2Q2(p + g)2 

simply by dimensional analysis. (There is no scale to make up for the fi"^"^.) So 



(F5) 



3 j: d'^'^p d'^^^q 



+ 



.\p + qI"^ - {p + qY \p - q\'^ - (p - 

+ (T independent) 



0, 



(F6) 
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where we have used the standard contour trick to do the sums and n{p) is the usual 
Bose function (|B3|) . The result is zero because (1) the temperature-independent piece in 
the penultimate line vanishes by dimensional analysis in 4 — 2e dimensions, and (2) the 
two terms 1/(P + Q)"^ and 1/(P — Q)^ in brackets exactly cancel each other. Because of 
this cancelation, the full result ( ^.391 ) for Isun is just equal to the purely three-dimensional 
contribution (p3D at leading order in the masses. 



2. Euclidean derivation 



In other computations in this paper, we will need to know the separate contributions of 
various subsets of Euclidean modes {po, go) to /sun- To get the formulas we need, we shall 
now rederive the result for Igun using purely Euclidean methods, similar to our derivation of 
the basketball integral. Start with 



^sun ~l~ -^sun ~l~ 



where Agnn and -Bsun are defined by 



rp2 



(477) 



In- 



mi 



mi + 1712 + ^3 



+ 0(m,e) 



(F7) 



A. 



2e 



^ n(P) 

> P2 ' 



^3 n(po=o,p) 

\3-2e 



(27r)^-^^ p^ + mf 

First consider Asun- As usual, we need to subtract out the UV divergences: 



(F8) 



(F9) 



A. 



(FIO) 



IP P2 TpP^TqQ'^ Tp 
where we have used the limiting behavior ( p. 27] ) of 11'^^) (P). The second term is now conver- 
gent in four dimensions. Exploiting the expression (p2|) for n(°)(P) and the integral (pT?!) 
enables us to evaluate the divergent parts of Asun as 



nW(P) 

P2 



{Any 



+ 2-7E 



T2 



(An) 



^ + iln^ 

12e 3 AttT 



0(e)) 

iC'(-i) 



3c(-i; 



'ATCfi"^ 

1 

+ 2 



+ 0(e) 



(Fll) 
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1 



T2 



1 2, /i IC'(-I) 1 1 



+ 0{e) 



(F12) 



fpP^TQQ^ (47r)2 

The finite part of A^^^ is calculated by utilizing the coordinate space integral representa- 
tion ( |2.25| ) for ]1^^\P) and ( p.30| ) for its UV behavior. Performing the p integral and the po 
sum produces 

T2 



(477)2 

(4^ 



roo \ / 1 r\ 

y dr- (cothr — 1) ( cothr ~ - ~ g ) ^(^) 

+ 0(e), (F13) 



2C' -1 1 , , , 1 
--^4 — f - -lE + In 27r - - 
3 C -1 3 ' ^ ^ 3 



where the method of appendix P have been used for the r integral. Adding ([Fll| , F12 , F13 ) 
yields 



A. 



^ n(P) 

p2 



T2 



(47r) 



— + In + ln(27r) + - 
4e 47rT ^ ^ 2 



+ 0(e) 



(F14) 



Now consider -Bsun- Though (|F9|) contains ultraviolet divergences due to the zero tem- 
perature part of n(P), dimensional analysis shows that the contribution from the zero 



temperature part is 0{m) and so need not concern us at leading order. Making use of (|2.25|) 
and completing the p integral, 



{An 



2 J A:7CT t'^ 



cothr — T ) + 0(m, e) 



(47r)2 Jo r 



( coth r — - — 1 I + 



(47r) 



r 



coth r — 



-mir 



+ 0(m,e) 
0(m, e) 



rp2 2^rp 

In h 0(m, e) 



(F15) 



(47r)2 mi 

where we have defined mi = mi/2nT and done the last step by the method of Appendix ^ 
Combining ^TS]) gives 



y2 



(4vr) 



1 , fL 1 

h In h - 

4e mi + 1712 + 2 



+ 0(m,e). 



(F16) 



Before leaving this section, we should collect some additional results that will be useful 
elsewhere. Subtracting (0) from ( P15|) gives 
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{P"^ + m?)Q^{P + Qf (47r)2 
Finally, adding ([1'12| ) and dFTS]) gives 



1 , 2T 1 

h In 

4e /i 2 



+ 0(m,e). 



(FIT) 



n(^)(p) _ 

p P2 - (4^^ 



6e ' 3 47rT 3 C(-l) 



1 + Hln " 



ln(27r) 



+ 0(e). 



(F18) 



APPENDIX G: EXAMPLE OF REDUCTION TO THE SCALAR BASKETBALL 



Consider the reduction of fig. |^(i): 



P-{Q-K){P-K)-Q 



8 "'"''fpQK p^qm^p - QY{Q - KY{K - py ' 

By expanding the numerator as in (|3.9| ), we get 



(Gl) 



{P-KYQ 



+ 



{P-KYQ 



+ 



P^Q^K^{P-Qy{Q-K)^ P^Q'^{P-Qf{Q-KY{K-PY 
[P-Kyq {P-K)-Q 



P^Q^K^{Q-Kf{K~Pf P'^K^{P-QY{Q-Ky{K-Pf 
Now switch the variables K and Q in the second lines: 



(G2) 



Id jpqk 



PiQ-K) 



PiQ-K) 



P^Q^K^{P-QY{Q-Kf P^Q^{P-Qf{Q-Ky{K~Pf_ 

(03) 



Use the identity ( |3.1(J| ) to substitute P, in the first numerator and K-^{P+Q) /2 in 

the second: 



64 JPQK 

^M4pqk 



1 



2P-Q-2P2 



P^K^{P-Qf{Q-Ky P^Q^{P-Qf{Q-Kf{K-Pf 
1 1 



P^K^{P-QY{Q-KY P^Q^{Q-Ky{K~Py 
1 



(G4) 



fpQK P^Q^K^{P + Q + K)^ ' 

where we have written 2P-Q as P^ — {P —QY for the second step and shifted integration 
variables in the last step. 



36 



APPENDIX H: DERIVATION OF I^^d 



The one-loop self-energy of fig. | can be reduced, using the methods discussed after ( p78| 
to the form 



where 



d-2 



qQ^{P + Q)^ 



(HI) 



n^,^ happens to be the form the self-energy would take in scalar QED. We find it convenient 
to introduce ll^j, mostly for reasons historic to our original derivation and because the 
decomposition simplifies some of the algebra of the following calculation. 

By again applying the same reduction methods, one may easily verify that both ([Hl|) and 
( P^ ) share the property that Pfjl^u = 0. In finite temperature non-Abelian gauge theory, 
this is a property of the one- loop self energy which does not persist to higher loops [|1^]. We 
shall use this property in our derivation. 

1. Consequences of Pj^ill^u = at one loop 

The orthogonality of Ufj^^ to P^ implies that it can be decomposed into separate transverse 
and longitudinal pieces:^^ 

n^,(p) = nT(P) Pt,.. + iIl{p) Plm. , (H3) 

where the Euclidean projection operators are given by 

VTij = Sij - PiPj/p^ , "Ptoo = VTOi = Vtm = , (H4) 
Plm. = S^u - P^^Pu/P^ - Ptm. . (H5) 
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For a review, see refs. 
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then gives 



(n 



+ (rf - 2)n2 



p2 
p2 



noo 



+ 



d-1 



' p2 ^ 



(H6) 



2. Scalar QED 



We now work to evaluate the integral 



-^sqed — 



p4 



and start by separating out the zero-temperature piece of Il^y 
T -4^ 

-fsqed — p4 



Anf;(p)]' + 2f -^nW(P)Ang:)(P) + ^ -1 [n(°)(P) 



(H7) 



(H8) 



a. The [finite temperature^ piece 

Let's evaluate the Po 7^ part of the sum for the first integral. First apply the standard 
reduction tricks to obtain 



n,,{P)=P'Tl{P) + 2{d-2)f^ — 



2 ' 



(H9) 



where n(P) is the scalar integral (p.l9|) . Next, we need to isolate the UV divergence of the 
P integration by isolating the large P behavior of H^,^: 



"S'p-2(d-i)|:^ 



1 



n 



(T) 



(HIO) 



00 P^oo P2 f Q2 

Specializing to the finite-temperature pieces of fl^j,, ( P6|) can be algebraically rewritten as 



2 _ {d-l) P' 

~ {d-2) p4 



n 



(T) 
00 



{d - 2) 



{d-2y 



n(^) 



(Hll) 



where we have introduced the UV subtracted 
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fi£'(P)^n£'(P)-^^i. 
n<T,(p)^n<-)(P)-A^i5. 



(H12) 
(H13) 



The integral we want is now 

^(T)Y _ {d - 1) ^' 1 r^(T) 



1 

P4 



The integrals in the last three terms can be obtained from ( p.l3| , p.33| , p.35| , |F18| ). We need 
to focus on the first two terms, which are convergent and may be evaluated with e = 0. 
From the observation that 



V 90 / 90 

where means (P /dr"^ and not V^, one may easily relate lloo to the scalar case 



1), (H15) 



T~tI {^-- [e''""''^(cothf + IpoD] - ie-|P»llPo|(p^ - 1)} 



+0{e) 

r 1 



e-l^'ol" fcothf- - - - 
r 3 



+ 0(e). 



(H16) 
(H17) 



Performing the angular integration and then integrating by parts yields 

The first two integrals in ( [H14| ) are then 
2 3 



2T p4 



n, 



(T) 
00 



2 (47r) 



POT^O 



oo poo 



JO 



J2 cise-l*'"l(^+'^M coth 



L r ) ( coth s 



r 3 



1 f d^p „r, sin nr sin ps 
— o: 



+ 0(e) 



(47r)2 



(27r)^ J p'^ ^ pr 

y Hdr rdse-^P°\^'+'^ ( cothf - - - - ) ( coths - 



1 r 



X ■ 



d^p ^2 sin P'" sin ps 
{2txY J p^ ^ pr ps 



0{e 



(H18) 

s 3/ 
(H19) 

s 3/ 
(H20) 
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Now plug in 



1 [ (fp „9 sin pr „9 sin ps 
/ o o 

(27r)3 J p4 



pr ps 
d^P_Q2 sii^ sin ps 



+ - — ^S{r-s) 



"(r - s) --^S{r - s) 



(H21) 
(H22) 



{2tt)^ J P^ ^ pr ps 27rr^ Arcr"^ 
Amazingly, the terms not proportional to 6{r — s) cancel in the sum of ( |H19| ) and ( |H20[ ). 
The 6{r — s) terms then give: 



2f p4 



n 



(T) 
00 



^'^n£)fi(-) = l^'[fi(-f + o(e). 



(H23) 



This could be easily evaluated using the techniques of appendix y, but we'll leave it in this 
form for now. 

The evaluation of the Po=0 piece of ( [FT^ ) proceeds in much the same way, but we don't 
need to make any UV subtractions. One finds 

^3 



Ang\0,p) = / d\ ^ (e'P-''- l) 92 (^cothf - ^) + 0(e) 

Tdrldf^] rcothf-iUo(e). 
JO \ pr J \ rj 



(H24) 



The integral is 

r I _fP_]_ 

(27r)3p4 



d^p f 3 1 r . - 



(27r)3 [2p^ 
+0(e) . 



An, 



(T) 
00 



2 1 



p 



(H25) 



The same sort of cancelation occurs between the first two terms as in the Po ^ case, and 
we are left with 



p4 



Po 



+ 0(6). 



(H26) 



Putting this together with the Po 7^ results (|H14|) and ( P23|) yields, after a little reorgani- 
zation, 

t j-. [^nf)] ' = f [n(^)] % 4(rf - 2)^ ^ ^) ' + 0(.) . (H27) 



The first integral is given by ( |2.35|) . One wonders if there's an easier way to get (P27|) 
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b. The rest of it 



The cross-term between n|^p and H^'J] is easy because U^^} is proportional to P'^S^y—P^Py. 
Using dHI), 



p4-^Vi^ -^Vi/ d _ p2 ^^^^ 



(H28) 



The integrals can be found in ( p. 13 , D2C , Ib'lll) . The final integral we need is 



p4 



n(o) 



d 



[n(o) 



(H29) 



which may be found in ( p2 Ij ) . Combining ( pi8| , H27 , [H28| , H29 ), and incorporating the 
results for the assorted basic integrals, gives 



1 



'rp2 \ 2 r 



'sqcd 



(47r)2 \12 



44C'(-3) , 272C'(-1) 



46 , u 

— + 92 In— . . , . . 

3e 47rT 3 C -3 3 C -1 



+ 167E 



1034 
1^ 



+ 0(e). 
(H30) 



3. Non-abelian gauge theory 



Using ([Hl|) and our standard reduction tricks, it is easy to obtain 

2 



d-2 



'qcd 



/,qed + 2rf/ball " A{d - ^n(P)^ ^ 



which, when combined with ( p.l3| , [F14| , |H30|) , is our final result for l^^A- 



(H31) 
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